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Abstract 

A study of the Instrumental background in balloon-borne gamma-ray 
spectrometers Is presented. The calculations are based on newly available 
interaction cross sections and new analytic techniques, and are the most 
detailed and accurate published to date. Results compare well with 
measurements made in the 20 keV to 10 MeV energy range by the Goddard Low 
Energy Gamma-ray Spectrometer (LEGS). The principal components of the 
continuum background in spectrometers with Ge detectors and thick active 
shields are 1) elastic neutron scattering of atmospheric neutrons on the Ge 
nuclei, 2) aperture flux of atmospheric and cosmic gamma rays, 3) e~ decays of 
unstable nuclides produced by nuclear interactions of atmospheric protons and 
neutrons with Ge nuclei, and 4) shield leakage of atmospheric gamma rays. The 
Improved understanding of these components leads to several recommended 
techniques for reducing the background. These include minimizing the passive 
material inside the shield and reducing the level of the shield threshold. A 
new type of coaxial n-type Ge detector with its outer contact segmented into 
horizontal rings can be used in various modes to reduce background in the 20 
keV to 1 MeV energy range. The resulting improvement in instrument 
sensitivity to spectral lines is a factor of * 2 in this energy range. 
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1. Introduction 

Due to the fntense radiation environment at balloon and satellite 
altitudes and the weak source strengths of astronomical objects at gamma-ray 
energies, spectropscoplc observations are necessarily background dominated in 
this energy range. Therefore, In addlton to increasing detector volume, one 
• can improve the sensitivity of gamma-ray spectrometers by decreasing the 

background level. One of the major challenges In designing new Instruments is 
to find the optimum detector and shield configuration and to make the right 
materials choices to achieve a minimum background. To facilitate such design 
efforts, we have performed a study of the sources of background in balloon- 
borne Ge spectrometers operating in the 20 keV to 10 MeV energy range. The 
accuracy of this work was improved over previous studies by using several new 
measurements of atmospheric gamma-ray fluxes and nuclear interaction cross 
sections. The level of success of the calculations was determined by 
comparing to background measurements made during two balloon flights from 
Palestine, Texas in 197!' and 1980 by the Low Energy Gamma-ray Spectrometer 
(LEGS) [1], an instrument developed at NASA/GSFC in collaboration with groups 
at CENS and Rice University. In this paper we present the results of this 
study, and describe several ways in which instrumental background can be 
reduced in future instruments. 

In general, the background in a balloon-borne gamma-ray spectrometer is 

ft 

made up of discrete background lines superimposed on a continuum. The origin 
t of the lines is natural radioactivity of the materials in the instrument, 

activation of the instrument materials by atmospheric neutrons, and the 
annihilation of positrons produced in the instrument and the atmosphere by 
atmospheric radiations. The origin of the continuum is atmospheric and cosmic 
gamma rays that enter the instrument through its aperture or penetrate its 
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shield, and activation of the instrument materials by atmospheric neutrons and 
protons. Most of the lines In the background spectrum are intrinsically 
narrow, and therefore appear in the measured spectra with widths equal to the 
instrument resolution. This ranges from several keV FWHM for high-resolution 
Ge spectrometers such as LEGS to several tens of keV FWHM for Instruments with 
scintillation-type detectors. 

Many examples of observed high-resolution background spectra and line 
Identifications can be found in the literature [1-7]. However, there are only 
a few published studies of the components that make up the background. For 
the continuum background there are three papers dealing with a phoswich 
scintillation detector [8-10] and one concerning a solid-state Hglg detector 
[113. For the 511-keV background line there Is one published study by Ling et 
al . [12] using data from a Ge detector. In this paper we present a detailed 
study of the continuum background In high-resolution Ge spectrometers and 
discuss techniques for its reduction. Preliminary results from this study 
have been given by Paciesas et al. [1] and Gehrels et al. [13], 

2. Instrumentation and Observations 

The LEGS instrument is described by Paciesas et al. [1]. In this section 
we will briefly review its general characteristics and give some details that 
are relevant for the background calculations. The Instrument performs high- 
resolution spectroscopy between -20 keV and 8 MeV using two interchangeable 
arrays of Ge detectors; the array chosen for a particular balloon flight 
depends on the observational goals. One array has three planar detectors of 1 
cm thickness, 53 cm 2 total effective area below 100 keV (as viewed through the 
aperture collimator discussed below), and 57 cm^ total active volume. Their 
efficiency is approximately unity between 20 and 100 keV, but drops steeply 
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above 100 keV. The other array has three coaxial detector of - 230 cm^ total 
active volume and 35.5 cin 2 peak effective area at 130 keV. The effective area 
decreases below 70 keV due to the top 1.4 mm Ge deadlayer typical of p-type 
coaxial detectors, and decreases above 200 keV due to the finite thickness of 
the detectors (4.6 cm). The decrease above 200 keV is less steep than that of 
the planar array, with the coaxial effective area at 1 MeV still 7.3 cm 2 . 

For both arrays, the Ge detectors are cooled to liquid nitrogen 
temperatures. This is accomplished by enclosing them In an evacuated cryostat 
and thermally connecting them to a dewar of liquid nitrogen by a copper 
coldflnger. The detector and cryostat are situated Inside an active Mai 
scintillator well to shield against cosmic rays and atmospheric radiations. 

The shield is - 13 cm thick, and collimates the field-of-view of the detectors 
to 16° FWHM. For the planar detectors, a passive Fe collimator is Inserted 
inside the shield to additionally collimate the field-of-view to 5° x 10° 

FWHM. The aperture solid angle per detector is given in Table 1 for the two 
LEGS flights for which background measurements are presented In this paper. 
Also listed In the table are the thresholds that were set on the shield pulse 
heights for the flights. Generally speaking, any photon or particle that 
deposits more energy than the threshold anywhere in the shield within - 2 ys 

A 

of a detector event vetoes that event. For the actual shield, nonuniformities 
in light collection from different regions of the Nal scintillator and 
suspected baseline shifts in the Nal photomultiplier tube amplifiers during 
the flight raised the effective thresholds in some regions by as much as a 
factor of w 2. 

Observations are performed by alternating every 20 minutes between 
pointing the detectors at a source and pointing at a region of sky without 
sources for background determination. The background measurements used for 
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comparison with calculations In this paper were obtained from background 
pointings in the time Intervals listed in Table 1. The ranges of atmospheric 
depth and pointing angle for the Instrument In these time Intervals are also 
listed In the table. 

As will be apparent in Section 3, an Important contributor to the 
background In present instruments is the passive material Inside the shield 
near the detectors. For LEGS, this Includes the cryostat, the detector 
housings, Inactive regions of Ge In the detectors, the section of coldfinger 
Inside the shield, the Inner walls of the shield housing, and the passive 
collimator for the planar detectors. Each material and its mass Is listed In 
Table 2. Another quantity used in the calculations Is the total Instrument 
material In the field of view. For the coaxial array in flight III, this 
Includes a 0.64 cm thick plastic scintillator at the top of the collimator and 
a 0.051 cm thick A1 window on the pressure vessel. For the planar array in 
flight V, there was no scintillator, and the A1 window on the pressure vessel 
was reduced to 0.038 cm. Both arrays also have material equivalent to 0.21 g 
cm -2 of air in their fields of view due to the air in the pressure vessel and 
the Insulation around the pressure vessel. The cryostat window thickness is 
0.17 cm A1 for the coaxial array and 500 pm Be for the planar array. The 
effective area plot in Paciesas et al . [1] includes the cryostat windows, but 
none of the other materials. 

3. Background Components 

In this section we discuss in detail the various components of the 
instrumental continuum background. The contribution of each component to the 
background of the LEGS instrument is calculated and compared to the measured 
background in flights III and V. Figures 1 and 2 show the measured background 
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spectra for the two flights, obtained with Instrument parameters and time 
Intervals listed In Table 1. The data Include only those events for which 
there was no simultaneous shield veto, since events with vetoes are excluded 
from source observations and therefore do not contribute to the background for 
the source measurements. As discussed in Section 1, a large number of 
background lines can be seen In both spectra. The anticoincidence shield 
surrounding the detectors suppresses the escape peaks and Compton continuum 
from these lines, although there is a small contribution as discussed In 
Section 3.5. In this paper we are calculating the level of the continuum 
background in the spectra, which is approximately equal to the background in 
the intervals between the lines in the measured spectra. The solid, dashed, 
and dotted curves in Figures 1 and 2 show the results of the calculations 
described below. 

3.1 Aperture Flux 

At low energies (< 150 keV), a major source of background in both the 
coaxial and planar arrays is atmospheric and cosmic gamma rays that enter the 
instrument aperture and Interact in one or more of the detectors. The photon 
need not deposit all of its energy In the detectors, but may also interact in 
passive material near the detectors, or even in the shield itself as long as 
the energy deposited does not exceed the shield threshold. 

The contribution from gamma-ray aperture flux was calculated using a 
modified version of the UCSD Monte Carlo photon/electron transport code [9] to 
simulate the interaction of gamma rays with the LEGS instrument. The incident 
downward gamma-ray flux over Palestine, Texas at 5 g cm” 2 for the coaxial data 
and 3.5 g cm" 2 for the planar data was estimated from measurements that have 
been made at balloon altitudes over the last decade. The available data at 5 
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g air 2 Is summarized In Figure 3. The spectrum used for the calculations Is 
Indicated by the solid line, and represents power-law fits to the combined 
data of Klnzer [14], Klnzer et al. [IS], and Schtfnfelder et al . [16], A 
similar analysis at 3.5 g cm -2 gives power law fits (ph cm" 2 sr" 1 s" 1 MeV" 1 ) 
of 2.19 x 10 2 E 0 * 70 between 0.024 and 0.035 MeV and 5.16 x 10" 2 E" 1 * 81 between 
0.035 and 10 MeV. The downward gamma-ray flux was used In this analysis, even 
though the pointing angle of the telescope during the background measurements 
reached values as large as 54° from zenith (see Table 1). However, this 
should result In only a small error since the total gamma-ray flux is 
approximately Independent of angle between 0° and 70° [18, 21]. 

The calculated contributions of the aperture flux to the continuum 
backgrounds for the coaxial and planar detectors are shown by the long-dashed 
curves in Figures 1 and 2 respectively. For the coaxial detectors there is an 
uncertainty of ±15% in this component due to uncertainties in our knowledge 
of the thickness of the inner Al wall of the shield housing (see Table 2) and 
of the actual shield threshold during the flight (100-200 keV). For the 
planar detectors, the same uncertainties in the shield housing and threshold 
exist, but the resulting uncertainty in the aperture flux curve from this 
source is only approximately ± 5%. This is because the Fe collimator tubes 
used in the planar configuration are the dominant passive material inside the 
shield and because the tubes tend to block low-energy secondary photons from 
interacting in the shield. 

The aperture flux curves in Figures 1 and 2 are broader and higher than 
one would calculate assuming a 0 keV shield threshold and no passive material 
inside the shield. This effect will be discussed in detail In Section 4.1. 

For the coaxial detectors the aperture flux is the largest component between 
70 and 150 keV, and for the planar detector It is the largest component 
between 30 and 100 keV. 
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3.2 Shield Leakage 

The other background component due to gamma rays Is the shield leakage. 

In this case, the background Is caused by the small fraction of atmospheric 
gamma rays that manage to leak through the shield without depositing more 
energy In It than the shield threshold, and then Interact In the detectors. 

The shield leakage contribution to the coaxial and planar detector background 
was calculated using the photon/electron transport code to simulate the 
Interactions of the Incident gamma rays with the Instrument. 

The zenith angle distribution of the total gamma-ray flux at MeV energies 
in the atmosphere has been shown by Schffnfelder et al. [21] and others to be 
essentially flat between 0° and 70° (0° = downward moving photons), rising by 
a factor of - 4 to a peak at 110°, and then falling by a factor of - 2 to a 
plateau between 130° and 180°. The actual spectra used as Input for our 
calculations were based on the measured downward fluxes described In Section 
3.1 and on measurements at 112° and 164° obtained at 2.5 g cm" 2 atmospheric 
depth by the Max Planck gamma-ray telescope [21]. The measurements at 2.5 g 
enr 2 were corrected to 5 g cm " 2 (coaxial detectors) and 3.5 g cm" 2 (planar 
detectors) with the calculated depth dependences of Ling [18]. Specifically, 
the forms for the differential spectrum (ph err 2 sr"* s"* MeV"l) used in 
various angular regions at 5 g cm -2 were 0.059 E"** 75 for 0° to 65° (Figure 
3), 0,094 E" 1 ' 61 for 65° to 95°, 0.15 E" 1 - 47 for 95° to 130°, and 0.047 
£-1.45 f 0r i 3 o° to 180°. At 3.5 g cm" 2 they were 0.052 E" 1 - 81 for 0° to 65*, 
0.085 E" 1 - 66 for 65° to 95°, 0.14 E" 1 - 50 for 95° to 130°, and 0.047 E“l- 48 for 
130° to 180°. 

There are several uncertainties In these spectra. The power-law forms 
were assumed in our analysis to apply over the energy range 0.1 to 10 MeV, 
while they were measured only In the 1.5 to 10 MeV range [21]. We feel this 
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Is not a major problem because, In the energy range of most uncertainty (below 
* 0.5 MeV), the dominant shield leakage background Is due to higher-energy 
photons that are scattered In passive materials near the detectors; the direct 
leakage of atmospheric photons Is cut off below 0.5 MeV by the Increasing 
absorption cross section of the shield. Another problem Is that a more recent 
flight of the Max Planck gamma-ray telescope [16] Indicated that the spectra 
presented In their 1977 paper (ref. 21; i.e., those used In our analysis) were 
low by a factor of 1.5 to 2 above - 5 MeV due to an overestimate of the 
background In that energy range. Unfortunately, angular distributions are not 
given for the later flight. It Is therefore possible that the shield leakage 
background component In the present calculations Is low above v 5 MeV. 

I»n the calculations, the simplifying assumption was made that the 
instrument pointed straight upward throughout the background measurements, 
although It was actually pointed more typically 20° to 40° from zenith. We 
expect the Inaccuracy in the results due to the pointing to be small (<< 30%), 
because there Is a cancellation of effects. As the Instrument tilts, the 
detectors see an Increase in shield leakage through one side of the shield as 
the thin upper section moves between them and the large flux at 90-110°, but 
at the same time they see a decrease through the opposite side of the shield 
as the thick lower section moves to 90°. 

The shield leakage component Is shown by the short-dashed curves In 
Figures 1 and 2. The uncertainty In the plotted curves Is estimated to be ± 
50% at 40 keV and ± 20% at 1 MeV for the coaxial detectors (Figure 1), and ± 
factor 2 at 40 keV and ± 40% at 1 MeV for the planar detectors (Figure 2). 

The shield leakage is seen to be a major component over most of the energy 
range of both detector arrays. For the coaxial array, the shield leakage Is 
virtually the only component above 2 MeV. The fit to the observed background 
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spectrum In this energy range Is seen In Figure 1 to be very good, Indicating 
that the shield leakage component Is well determined by the above calcuation 
technique. The agreement with the data Is In fact fortuitously good given the 
Input spectrum uncertainties and uncertainties In the shield threshold and 
wall thickness. As In the case of the aperture flux, the nonzero shield 
threshold and the passive material Inside the shield cause an Increase In the 
background, as will be discussed In Section 4.1, 

3.3 Beta Decay 

The 3-decay background Is due to the decay of 3"-unstable nuclides 
produced by nuclear Interactions of atmospheric protons and neutrons with Ge 
atoms in the detectors. The signal In this case is produced by the ionization 
energy loss of the decay electron as It stops In the detector. Since the 
typical °“-decay lifetimes are much longer than the microsecond coincidence 
•vi.ms-s ‘if the Instrument, any interactions of the primary proton or neutron In 
the shield do not veto the event, and any energy deposition in the detector 
from the primary particle or the recoil of the Ge atoms is not Included in the 
signal. 

In many cases, the daughter nuclide of the 3 decay is produced in an 
excited state that decays in a time short compared with the coincidence 
time. Therefore, In addition to an electron and neutrino, these decays also 
have one or more prompt gamma rays and occasionally (Internal conversion) 
another electron. For the case where no gamma rays are emitted (i.e., decays 
to the ground state or a metastable state of the daughter), the energy is 
deposited in a small localized region of the detector due to the small range 
of the decay electrons; for example, the range of a 1 MeV electron in Ge Is 




< 1 mm*. For the case where one or more prompt gamma rays are emitted, the 
event will fall in one of the following categories: 1) The gamma ray does not 

interact in the detector and does not deposit enough energy in the shield to 
exceed the threshold. The event is therefore accepted, and appears as a 
localized event. 2) The gamma ray deposits enough energy in the shield to 
exceed the threshold, causing the event to be vetoed. 3) The gamma ray 
interacts in the detector, but does not deposit enough energy in the shield to 
exceed the threshold. The event is accepted and the signal in the detector Is 
the sum of the gamma-ray deposition and the (T energy loss. Since the gamma- 
ray interaction Is, in general, in a different region of the detector than the 
3" decay, these events are called nonlocalized events. 

Only 3~ decays are discussed in this section, although 3 + decays can also 
produce background events. The ccritinuum background contribution from 3 + 
decays in the detector was calculated and found to be extremely small, 
primarily because the events are suppressed by the photons produced when the 
positron annihilates. The annihilation produces two oppositely-directed 511- 
keV photons, at least one of which escapes the detector in almost all cases 
and Interacts in the shield vetoing the event. The decays that proceed via 
electron capture do not produce positrons, and are therefore not vetoed by the 
annihilation radiation. However, these decays generally do not contribute to 
the continuum background since any prompt gamma rays absorbed in the detector 
will be In narrow background lines. 

Since the 3~-deca.y lifetimes a rc much longer than che instrument 
coincidence times, the relevant particle spectra at the detectors should 
include all contributions from secondaries produced in the shield. In order 
to obtain these spectra, we followed the suggestion of Mahoney et al. [6] of 

*The electron ranges used in this paper were calculated with formulas in Evans 
[22] for extrapolated ranges in Al, scaled to other materials using the 
relation R(A,Z) = A/Z *13/27 *R^-j given by Koral and Cohen [23]. 
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using spectra measured In the atmosphere at a depth equal to that of the 
actual atmosphere above the instrument plus the equivalent depth of the 
shield. The LEGS shield is approximately 48 g cm" 2 thick, which corresponds 
to 31 g cm -2 of air due to the difference in nuclear collision lengths (60.2 g 
cm" 2 for air vs 91.1 g cm" 2 for Nal [24]). Thus, with actual instrument 
atmospheric depths for the LEGS flights in the range 3 to 5 g cm" 2 , the total 
effective atmospheric depth inside the shield is ^ 35 g cm" 2 . 

The Input neutron spectrum for this calculation was derived from those 
calculated by Armstrong et al. [25], which are consistent with observations 
and are presented in a convenient format. The specific spectrum we used is 
shown in Figure 4, and is a logarithmic interpolation to 35 g cm" 2 between 
power-law representations of their spectra at 10 and 50 g cm" 2 . The proton 
spectrum below 3.6 GeV (geomagnetic cutoff at Palestine) was more difficult to 
obtain because there are no recent calculations or published measurements. 

The approach taken was to use the equations of Rossi [26] for the atmospheric 
proton and neutron spectra, and determine the intensity constant in the 
equations by fitting the neutron spectrum to the one shown in Figure 4. The 
resulting proton and neutron curves are shown by dashed lines in the figure. 
The shape of the Rossi neutron spectrum is seen to be in reasonable agreement 
with that of the Armstrong et al. spectrum. The proton spectrum used in the 
3" calculation is shown by the solid-line power-law curves, and was derived by 
matching the Armstrong et al. neutron spectrum above 300 MeV and following the 
shape of the Rossi proton spectrum at lower energies. There is one set of 
measurements by McDonald and Webber [27] of the integral proton flux between 
100 and 750 MeV at a similar magnetic latitude as that of Palestine. The flux 
at 35 g cm -2 implied by their observations is 0.23 protons cm" 2 s"l assuming 
an isotropic flux distribution, which is in excellent agreement with the value 
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calculated from the power-law curves in Figure 4 of 0.25 protons enr*- s" 1 . 
Above the geomagnetic cutoff at 3.6 GeV, the proton spectrum was taken to be 
the primary cosmic-ray spectrum measured by Webber and Lezniak [28] multiplied 
by a factor of 0.56 to account for protons lost to nuclear collisions in the 
atmosphere and in the shield. 

An extensive literature search was performed to find the required 
interaction cross sections for protons and neutrons incident on the five 
naturally occurring Ge isotopes to produce ^“-unstable nuclides. For 
neutrons, several compilations of cross section data are available [29-32]. 

By far the most useful of these for our application was the compilation of 
Howerton, Dye, and Perkins [32] which has recently been revised to include 
data on Ge. For protons, the semi empirical cross sections of Silberberg and 
Tsao [33,34] were used. The calculations were greatly facilitated by a 
computer code provided to us by the authors giving their most recent cross 
sections. The accuracy of the present study would have been significantly 
poorer without the new neutron cross sections for Ge of Howerton, Dye, and 
Perkins, and without convenient access to the proton cross sections of 
Silberberg and Tsao. 

The actual calculations are described in detail in the Appendix, and will 
only be summarised here. The ^'-activation rates were determined by 
integrating the neutron and proton spectra at the detectors with the 
interaction cross sections of each Ge isotope. All final states listed as B” 
emitters In the Table of Isotopes [35] were considered, amounting to more than 
a hundred for each Ge Isotope. Each activation rate was then multiplied by 
the 3" branching fraction and the decay mode fraction for each possible final 
state, to give the rates for the different dbv odes of each nuclide 


produced. Decays accompanied by a prompt gamma ray were also multiplied by 
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another factor to eliminate events with Interactions In the shield. It was 
assumed that any gamma ray not fully absorbed In the detector Is vetoed by the 
shield. The efficiencies for intenal ly-produced photons to be fully absorbed 
In the detectors were calculated using the Monte Carlo code. For LEGS this 
assumption gives an underestimate of the background produced by S" decays 
accompanied by prompt gamma rays, since some of the photons escaping the 
detectors are absorbed In the passive materials near the detectors or Interact 
In the shield without exceeding the high (~ 100 keV) threshold. However, tbe 
effect on the total background is negligible, because, as will be shown, the 
prompt gamma-ray decays are only a small fraction of the background. In later 
sections of the paper, these calculations will be repeated for future 
Instrument configurations that have very little passive material near the 
detectors and low thresholds for the shield. In these cases the use of the 
total-absorption efficiency for the nonvetoed fraction Is nearly correct, and 
the calculated rates for the 3“ decays accompanied by prompt gamma rays should 
be accurate. 

For each decay mode for each ^“-unstable nuclide produced, the 3" energy 
spectrum was determined using the rates calculated as described above and the 
3" spectral shapes of Behrens and Szybisz [36]. The spectra were then summed 
to produce the total 6"-decay background spectrum in the detectors. More than 
1000 production rates were calculated and 60 spectra summed to determine the 
final background spectrum. The results for the LEGS coaxial and planar 
detector arrays are shown by the dotted curves In Figures 1 and 2 
respectively. For the coaxial array the background produced by the 
nonlocalized prompt gamma-ray decays is shown separately from the localized 
ground-state decays. For the planar array only the localized decays are 
Included, because the prompt gamma rays from the nonlocalized decays typically 
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escape the thin detectors and interact In the shield, vetoing the event. 

One of the new results obtained from this calculation is that the 
background from nonlocal ized decays is extremely small - more than an order of 
magnitude less than the background from the localized decays. This Is a 
somewhat surprising result given that ~ 2/3 of all 0~-decay nuclides that can 
be produced by neutron and proton interactions with Ge have a greater than 50% 
probability of emitting prompt gamma rays during the decay. The reason for 
this is twofold: first, by chance, the three most abundantly produced 

nuclides, namely 75 Ge, 70 Ga, and 69 Zn {see Table 4 and Figure 18), decay 
predominantly to the ground states of the daughter nuclides; second, the 
decays with prompt gamma rays are suppressed because many of the photons 
interact In the shield or produce secondaries that interact In the shield and 
veto the event. 

The P“-decay component is an important part of the total background 
between - 150 keV and - 1,5 MeV for the LEGS coaxial detectors, and above 
300 keV for the planar detectors. For future instruments with thicker 
shields and smaller apertures than LEGS, the 3"-decay component will dominate 
the background at intermediate energies, as will be shown in Section 4.3. 

3.4 Elastic Neutron Scattering 

The elastic neutron scattering background is due to the recoil of Ge 
atoms in the detectors following elastic scatterings with Incident atmospheric 
neutrons. The signal from this interaction is prompt, so that only those 
neutrons that do not interact in the shield while entering or leaving the 
instrument contribute to the background. The calculation of this component 
was done by determining the spectrum at the detectors of neutrons that do not 
interact in the shield, and integrating the spectrum times the Ge elastic 
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scattering cross sections. The energy deposited In the detector was found by 
multiplying the recoil energy by the fraction of energy lost via signal- 
producing ionizing collisions in the detector. 

For the Incident neutron spectrum, we used power-law approximations of 
the calculated spectrum of Armstrong et al . [25] at 5 g cm" 2 atmospheric 
depth. This spectrum was multiplied by a correction factor to eliminate 
neutrons that veto the background event upon entering or leaving the 
instrument by Interacting In the Nal shield and depositing more energy than 
the shield threshold. To obtain this factor It was necessary to determine 
what kinds of interactions in the shield are capable of depositing more energy 
than the threshold. In particular, for elastic scattering, the recoil energy 
of the nucleus, E R , is given by 

e r“TW (1 - cos0) E n (1) 

(see, e.g., ref. 37) where A is the mass of the target nucleus divided by the 
neutron mass, E n Is the incoming neutron kinetic energy in the lab frame, and 
0 is the scattering angle of the neutron in the center-of-mass frame. For Na, 
the maximum possible recoil energy Is only 16% of the incident neutron energy, 
and, for I, only 3%. Thus, based on equation (1) alone, it can be seen that 
neutrons with energies less than 0.6 MeV can not exceed the - 100 keV shield 
threshold via elastic scatterings. 

It can be shown, in fact, that even neutrons with energies in the 1-100 
MeV range generally will not trigger the shield via elastic scatterings. This 
is because the elastic scattering cross sections at high energies are peaked 
near cos 0 = 1 [38] and because, as will be discussed below, in this energy 
range more than half of the recoil energy is lost to interactions in the Nal 
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that do not produce detectable signals. Inelastic collisions, on the other 
hand, generally disrupt the target nuclei producing photons and/or particles 
with typical energies in the MeV range, and therefore almost always produce 
shield vetoes. Hence, we make the assumption that no elastic scatterings 
exceed the shield threshold but that all inelastic scatterings do. The 
Interaction length due to inelastic scatterings (also know as the absorption 
length) In Hal Is 41.3 cm [24], so that the fraction of incident neutrons that 
penetrate the -• 13 cm thick LEGS shield twice (entering and leaving) without 
producing a veto is 0.53. In later sections of this paper, a new instrument 
will be considered with a 15.2 cm (6") thick Nal shield. For this case, the 
noninteracting fraction is 0,48. The actual neutron spectrum (neutrons cm -2 
s" 1 MeV -1 ) at the detectors used in the LEGS calculations, Including the 0.53 
factor., is 0.053 E n "'® 3 between 0.1 and 100 MeV and 0.85 E n " 1,6 between 100 
and 10 4 MeV. 

The equation for calculating the counting rate per unit volume of Ge, per 

dR 

unit energy detected signal in the detector, -gp for the elastic neutron 
scattering background Is 


dR _ dE R dR 
dE~W~ BE^ (E R" E/f) 


1 

f + df/dE R • E 


•R 


10" 27 p N, 


w 


/°°da , 

'oV" 


dE„ ents s” 1 MeV -1 cm -3 (2) 


where E Is the energy collected as signal in the detector, E R is the Ge atom 
recoil energy, f Is the fraction of the recoil energy collected as signal 
(= E/E r ), p is the density of the target element (5.36 g cm" 2 for Ge), is 
6.02 x 10 23 atoms mole" 1 , w is the atomic weight of the target element (72.6 g 
mole" 1 for Ge), is the elastic scattering cross section per unit recoil 
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energy for incident neutrons of energy E n in units of mUlibarns MeV"*, 10" 2 ^ 
Is the number of cm 2 per minibarn, and j n is the neutron flux at the detector 
at energy E n In units of cm -2 s" 1 MeV"* (see above discussion). The cross- 
section data used In this analysis are given in terms of as a function of 
o, so we take one more step, using equation (1) to convert from to 
with the result for Ge 


dR _ 1.03 x 10" 2 { m da 


(3) 


where the lower limit on the integral is the minimum energy neutron that can 
produce a Ge recoil of energy E R , determined from equation (1). 

The cross sections for elastic neutron scattering on Ge were obtained 
from the plots of Garber et al . [38]. The data are given in that compilation 
as plots of vs 0 for various E n . in order to convert to ^ vs E n for 
constant E R (or E), one point was taken from each E n plot at a value of 0 
determined using equation (1). The results for two example energies, E R = .04 
and .17 MeV, are shown by the filled circles in Figure 5. In both cases, 
there is not enough data on Ge to adequately define the shapes of the curves 
at energies greater than ~ 10 MeV. For this reason, we have used the more 
extensive data on Cu to indicate the shapes at higher energies. For each 
plot, the values of E R for Cu were chosen such as to have the same E n m ^ n as 
Ge. Specifically, we chose E R (Cu) = 1.138 E R ( Ge) (equation (1) with A(Ge) = 

71.97 and A(Cu) = 63.00). The cross sections for the Cu are shown as x's in 

r 

Figure 5. For small values of E R (Ge) near .04 MeV, the agreement between the 
Ge and Cu points is quite good in the energy range where both are measured, 
whereas for the larger values of E R (Ge) near .17 MeV, there are some 
systematic disagreements in the trend of the points for the two elements of 



typical magnitude 50$. The lines in Figure 5 are the power-law approximations 
to the cross-section curves that were used in the integral in equation (3). 
Approximately half of the contribution to the Integral is from E n > 1° MeV. 

The final consideration in the elastic neutron scattering calculation is 
the fraction of the recoil energy that produces signal in the detector. 
Qualitatively, an elastic scattering interaction in a Ge crystal can be viewed 
as follows: The neutron enters the detector and scatters elastically off the 
nucleus of a Ge atom imparting a recoil kinetic energy to the atom nearly 
equal to the energy lost by the neutron (the displacement energy of a Ge atom 
from the lattice is only ~ 18 eV [39]). The Ge atom then loses its energy 
both via electronic collisions that produce ionization in the detector 
(collected as signal), and via atomic collisions that displace other atoms 
from the lattice. The atomic collisions do not produce Ionization in the 
detector, although the displaced atoms can, in turn, lose part of their energy 
to electronic collisions. The total fraction of the recoil energy that is 
lost to signal-producing electronic collisions is denoted by f in equations 
(2) and (3). 

The fraction f has been determined theoretically by Lindhard et al . [40], 
with the relevant equations given in a convenient format by Robinson [41]. 

The equation used for the present analysis is 


L 9< E R/ E L> 


where k L = 0.134 Z 2 / 3 A" 1 / 2 , E L = 86.93 Z 7/3 eV, g(e) = 3.401c 1 / 6 + 0.402 
E 3/4 + e, z is the atomic number, and A is the atomic mass. Figure 6 shows f 
as a function of recoil energy, E R , for Ge. As an example using these data, 
we see that for the recoil energies of .04 and .17 MeV discussed in regard to 
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Figure 5, the measured energy in the detector is .012 and .066 MeV 
respectively. 

Combining the neutron spectra, the elastic scattering cross sections, and 
the data In Figure 6 with equation (3), gives the elastic neutron scattering 
background components in Figures 1 and 2 for the LEGS coaxial and planar 
detector arrays. The main uncertainties in these curves comes from 
uncertainties in the cross-section data such as those shown in Figure 5. We 
estimate this uncertainty to be 20 % at low energies {E R “ *04 MeV, E = .012 
MeV) and 70 % at higher energies (E R > .17 MeV, E > .066 MeV). The elastic 
scattering component is seen to be important at low energies, exceeding even 
the aperture flux below 50 keV for the coaxial detectors and below 25 keV for 
the planar detectors. 

3.5 Other Components 

We now estimate the magnitude of several other components and show that 
their contributions to the background are small. 

1. Continuum Due to Spectral Lines - The thick anticoincidence shield 
that is used to keep out atmospheric radiations in gamma-ray spectrometers 
also serves as an excellent veto for events with Compton-escape photons. 
However, the nonzero shield threshold and the passive material between the 
detectors and the shield, allow some such events to be included in the 
background. For the background components discussed so far, this effect 
has been included in the calculation, but there is an additional component 
due to the Compton continuum of the lines in the spectrum. We consider 
for example the line at 198 keV, which is the strongest line in the 
coaxial spectrum. 

The interaction producing the 198 keV line is 7°Ge(n,Y)7* m Ge. The 
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metastabl e state decays with a half-life of 22 ms yielding a 23 keV 
Internal conversion electron and a 175 keV photon. Using laboratory data 
and Monte Carlo simulations we estimate (to within a factor of 2) that the 
continuum background below a line of this energy is 4 cnts MeV" 1 per 
count In the line. There Is some structure In this continuum such as a 
Compton edge and a broad backscatter peak, which we ignore for the purpose 
of this approximate calculation. For the coaxial detectors In flight III, 
the intensity In the 198 keV line was 0.65 cnts s" 1 , giving a continuum 
background at energies less than 198 keV of - 1.1 x 10~ 2 cnts s" 1 MeV" 1 
cm" 3 . Comparing wfth Figure 1, it Is seen that the contribution Is 10-20% 
of the observed background in this energy range. This component Is not 
Included In Figure 1 because the data were not available to perform an 
accurate calculation. In future Instruments with thicker shields, less 
passive material near the detectors, and lower shield thresholds, this 
component should be negligible. 

2. Electron Aperture Flux - The electron aperture flux background Is due 

to primary and atmospheric-secondary electrons that enter the aperture of 
the Instrument and stop in a detector. We obtain the Incident flux In the 

upper atmosphere from the calculations of Daniel and Stephens [42], who 

treat both electrons and gamma rays with energies between 1 MeV and 10 GeV 

at all atmospheric depths. There Is some disagreement between their 
gamma-ray spectra and observations below 1 GeV [43], but their electron 

spectra agree well with observations [44]. At 5 g cm" 2 atmospheric depth, 
the calculated electron Intensity can be represented by 1.4 x 10“ 2 E" 1 * * 8 
electrons cm" 2 sr" 1 s" 1 MeV" 1 between 1 and 10 MeV, and is approximately 
Isotropic below 5 MeV [42]. 

To set an extreme upper limit on the electron contribution to the 
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background, we assume that the spectrum between 1 and 10 MeV can be 
extrapolated as a power law to lov/er energies, and that there Is no 
absorbing Instrument material In the aperture. If these assumptions were 
valid, the electron background would be significant at low energies; for 
example, with an aperture solid angle of .08 sr and a detector area of - 
50 cm^ the background for the coaxial detectors would be .04 cnts s" 1 MeV 1 
cm" 3 at 60 keV, which Is comparable to the aperture flux, although It would 
be only 2 x 10"4 cnts s" 1 MeV 1 ' cm" 3 at 1 MeV, which Is small. The above 
assumptions, however, are not valid. The spectrum almost certainly falls 
below the power-law extrapolation at energies less than 1 MeV, because the 
electron range In the atmosphere becomes small at these energies. Also, 
the material In the path of the Incoming electrons Is thick to low energy 
electrons. The 47 cm thick A1 window on the cryostat, by Itself, stops 
electrons with energies less than - 1 MeV. Including the A1 pressure- 
vessel window, the plastic scintillator In the aperture, and the inactive 
layer of Ge on the detectors, the threshold energy for electrons to reach 
the active Ge for the coaxial detectors Is over 3 MeV, Hence, the 
electron aperture flux Is a negligible background component. 

3. Proton Aperture Flux - The proton aperture flux background, due to 
atmospheric protons that enter through the instrument aperture and stop in 
the detectors, can be seen Immediately to be a negligible component,. 

First, the flux of protons Is extremely small; even at 35 g cm" 2 
atmospheric depth, the intensity at 1 MeV for an Isotropic flux Is « 10"° 
cm" 2 sr" 1 s" 1 MeV" 1 (see Figure 4) compared with 1.4 x 10" 2 for the 
electrons. Also, the proton range is less than that of electrons, so 
that, for Instance, even the .17 cm thick cryostat window will stop all 
protons with incident energy less than 18 MeV [45], 
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4. Background Reduction Techniques 

The sums of ... calculated components shown by the solid lines In Figures 
1 and 2, are seen to be generally consistent In both shape and magnitude with 
the observed continuum background spectra from the LEGS coaxial and planar 
detector arrays. Since the two detector configurations are quite different In 
terms of their fields of view, passive materials near the detectors, and 
detector volumes, the success of the calculations for both cases Indicates 
that the background components are well understood, in this section, we study 
In more depth the nature of the different components, and suggest techniques 
for their reduction. 

4.1 Passive Material and Shield Threshold 

The aperture flux and shield leakage curves In Figures 1 and 2 are 
significantly higher and broader than one would calculate for a very low 
shield threshold and no passive material Inside the shield. The extra 
background Is due to photons that either scatter in the detector and rtre 
absorbed in the passive material or In the shield { AE < shield threshold), or 
vice versa. 

For the aperture flux, this increase Is Illustrated In Figure 7, where 
the ratio of the total aperture flux background to the contribution from 
unscattered photons Is shown as a function of energy for the coaxial and 
planar detectors. The unscattered contribution Is that due to photons that 
enter the aperture and are fully absorbed In the detectors. It can be simply 
calculated at a given energy by multiplying the downward gamma-ray flux 
{Figure 3) times the detector effective area (Figure 2 of Paclesas et al . Cl]) 
times the aperture solid angle (Table 1) times an absorption factor 
(approximately equal to 1) to account for the instrument materials in the 
field of view (Section 2). 
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Two curves are shown In Figure 7 for the coaxial array, representing the 
extremes of uncertainty In the thickness of the Inside walls of the shield 
(0.16 to 0.32 cm A1 ) and the shield threshold during the flight (100 to 200 
keV). The coaxial curves are high at low energies because the unscattered 
spectrum falls off steeply toward lower energies due to the Ge deadlayer on 
the detectors, while the contribution from scattered photons Is relatively 
constant In this energy range. The ratio reaches a minimum near 70 keV and 
then rises toward higher energies because of the Increasing cross section for 
Compton scattering relative to photoelectric absorption. 

For the planar array, the ratio In Figure 7 Is approximately constant 
between 40 and 80 keV because the detectors have essentially no top dead- 
layer. The steep rise above 100 keV Is caused by the Increasing Compton- 
scattering cross section In this energy range In the Iron collimator. The 
factor of ^ 5 Increase In aperture flux background above 200 keV due to the 
massive Fe collimators close to the detectors in the planar array illustrates 
how significantly passive material Inside the shield can enhance the 
background. The collimators are still a valuable feature of the LEGS planar 
configuration, however, since they restrict the field of view and thereby 
decrease the background In the Important energy range below 100 keV. 

Table 3 itemizes the contribution of different passive materials Inside 
the LEGS shield to the Increase in aperture flux background In the 0.1 to 1.0 
MeV range for the coaxial array. Out of the 71* total Increase (thick shield 
housing, 200 keV shield threshold), 39* Is due to single scatterings in a 
passive material, 13* Is due to multiple scatterings, and 19* Is due to events 
with energy (< 200 keV) deposited in the shield. We have not found a simple 
quantitative relationship between the increase In background caused by a 
particular material and the material parameters, such as mass, atomic number, 
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thickness, distance from detector, and covering factor. However, It Is clear 
considering background Increase per unit mass, that passive materials very 
near the detector, such as the Ge deadlayer and the A1 detector housing, are 
particularly bad. The higher atomic number of the Ge also plays a role since 
the Ge deadlayer contribution Is 2 1/2 times that of the A1 detector housing, 
while the materials have similar masses and are both close to the detector. 

Now turning to the shield leakage background, the Increase due to the 
high shield threshold and passive material near the detectors is shown for the 
coaxial array In Figure 8. The curves labeled 1 and 3 are ^he total 
calculated shield leakage backgrounds for the two extreme cases of 0.32 cm 
thick shield housing walls and 200 keV shield threshold, and 0.16 cm thick 
walls and 100 keV threshold. The shield leakage curve used In Figure 1 is the 
geometric mean of these two curves, and is shown In Figure 8 by the dashed 
line, curve 2. The curve labeled "Unscattered Photons" represents the shield 
leakage background that would be obtained for an ideal instrument with only 
active Ge Inside the shield and a 0 keV shield threshold. The most dramatic 
effect the passive material and shield threshold has on the background occurs 
at energies less than ~ 0.4 MeV. The unscattered component falls off steeply 
due to the rising absorption cross section of the shield, whereas the total 
background remains high. In this energy range, the background is caused 
almost entirely by higher energy photons that either scatter in the shield 
without triggering the threshold or scatter in passive material near the 
detectors. Curve 4 in Figure 8 shows the effect of the passive material by 
Itself with a 0 keV shield threshold. At 1 MeV, the passive material inside 
the shield doubles the shield leakage background, and the 100 keV shield 
threshold contributes an additional 60%. 

In order to further explore the effect of shield threshold on the shield 
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leakage background, we have performed Monte Carlo simulations for an 
instrument configuration that may be typical of new instruments. The 
instrument was assumed to have a 15.2 cm thick Nal shield surrounding (except 
for a 20° FWHM aperture hole) several coaxial n-type Ge detectors of dimension 
7 cm diameter x 7 cm length**. No passive material was included inside the 
shield so that the effect of the shield threshold could be determined 
independently of other parameters. The calculation is for a flight at 3.5 g 
cm" 2 atmospheric depth over Palestine, Texas. The results are shown in Figure 
9, where the ratio of the total shield leakage background to the leakage with 
a 0 keV threshold is plotted as a function of shield threshold for three 
energy bands. At energies greater than 2.5 MeV, where the shield leakage is 
by far the dominant background component, the increase in leakage due to the 
nonzero threshold is only in the few percent range for reasonable thresholds, 
but does approach 50% for high ( v 200 keV) thresholds. At lower energies the 
increase is much greater, but the total effect is not as important because the 
3“ component is also a significant fraction of the background. However, as 
will be discussed in Section 4.3, reduction techniques can be used in new 
instruments to substantially reducing the 3" component. In this case, the 
increase in the shield leakage background at lower energies can lead to a 
significant increase in the total background. Based on Figure 9, a shield 
threshold well below 50 keV is recommended. 

The standard techniques for reducing the aperture flux and shield leakage 
backgrounds involve changing instrument parameters, such as aperture size and 
shield thickness. We have shown here that the background can also be 

**For all calculations in this paper concerning coaxial detectors other than 
the LEGS coaxials, it is assumed that the cylindrical hole in the center of 
the detector for the inner contact is 1 cm in diameter and stops 1 cm below 
the top of the detector. 
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significantly reduced without such design changes by decreasing the amount of 
passive material Inside the shield and by lowering the shield threshold as 
much as possible. At low and high energies where the aperture flux and shield 
leakage respectively dominate the background, a reduction In total background 
by more than a factor of 2 can be achieved, which means an Improvement in 
sensitivity to spectral lines of more than 40% (see equation (6)). 

Specific recommendations for reducing passive material and lowering the 
shield threshold are: 

1) Use n-type coaxial Ge detectors with the reverse electrode configuration. 
For these detectors the thick (1-2 mm) Li-diffused deadlayer of Ge can be on 
the inner contact rather than the outer contact as Is the case with the 
standard p-type detectors [46]. This reduces the volume of inactive Ge by a 
factor of - 9 assuming a coaxial detector of 7 cm diameter and 7 cm height. 

2) Reduce the mass of the detector housing and cryostat to an absolute 
minimum, and use low-Z materials such as A1 or. Ideally, Be. 

3) Replace the standard Cu coldfinger with an A1 coldfinger, which gives the 
same thermal conductance for approximately half the mass. 

4) Reduce passive materials In the aperture to an absolute minimum. This can 
be done by using an active collimator whenever possible. If a small field-of- 
view requirement dictates a passive collimator, the background produced by the 
collimator can be minimized by first using an active collimator, such as an 
aperture hole in the shield, to reduce the field of view to a few tens of 
degrees, and then placing the fine collimator at the end of the hole. The 
geometry factor for scattering Into the detector Is thereby much reduced 
compared with placing the passive collimator near the detectors. The inner 
housing walls of the aperture hole in the sheild also contribute to the 
scattering background. This contribution can be eliminated by using thin Be 
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windows on either side cf the hole, so that the inner housing wall is 
eliminated altogether. The Be will be directly in the field of view, but the 
attenuation of the incoming photon beam is extremely small; for instance, if 
windows of 0.05 cm thickness are used, the attenuation at 20 keV for two 
windows is only 4%. 

5) Use Na I scintillator for the shield in order to set as low a shield 
threshold as possible. Experience has shown that setting a low threshold, 
below 50 keV, in a large-volume Nal shield at balloon altitudes, without 
producing an unacceptable deadtime, is an extremely difficult challenge. Of 
the three commonly used scintillators (Nal, Csl, BGO), Nal has the highest 
light output per unit energy loss in the shield (relative numbers 1.0, 0.85, 
0.13) and the shortest light decay constant (0.23 us, 0.63 us, 0.3 us). It is 
therefore even more difficult to set low thresholds in the other two 
scintillators than in Nal. There are many other important criteria that enter 
into the decision of what scintillator to use for a shield (such as cost, 
availability, Instrument geometry, and weight), but strictly from the 
standpoint of minimizing shield leakage background by setting as low a shield 
threshold as possible, Nal is the scintillator of choice. Experience has also 
shown that light collection uniformity from the scintillator is a critical 
parameter in actually achieving a low threshold during a flight. 

All of the above recommendations are being incorporated into a new 
instrument that is currently being built as a collaboration among groups at 
Bell Laboratories, Goddard Space Flight Center, and Sandia National 
Laboratories, Albuquerque. The total mass of passive material inside the 
shield for this instrument is ~ 3.5 kg (almost all of which is Al), which is a 
similar mass to the ~ 5 kg for the LEGS coaxial array. However, the new 
instrument is designed for seven Ge detectors, each 7 cm in diameter by 7 cm 
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in length. Therefore, the total passive mass per unit volume of Ge Is - 2 g 
cm" 3 which Is an order of magnitude less than the - 20 g cm -3 for LEGS. 

In addition to decreasing the continuum background, reducing the passive 
material Inside the shield also results In a decrease In the Intensity of the 
background line at 511 keV. Part of this background is due to positrons which 
are created in the passive material via pair-production and nuclear 
interactions, and then annihilate producing two 511-keV photons. This energy 
is of particular interest for present and future gamma-ray spectrometers since 
positrons can also be produced In astrophyslcal sources. 

4.2 Aperture Size 

In the energy range from 30 to 100 keV where the aperture flux dominates 
the background, the most straightforward approach to reducing the background 
is to decrease the field of view of the Instrument. The effect can be seen by 
comparing the aperture flux for the LEGS wide fiel d-of-view coaxial array with 
the narrow field-of-view planar array. Since the relevant quantity to compare 
for front-incident low-energy photons is the background per unit detector 
effective area, we normalize the backgrounds in Figures 1 and 2 by the coaxial 
and planar^effectlve areas in Figure 2 of Paciesas et al. [1], At 100 keV, 
the coaxial aperture flux background per unit detector effective area is 0.33 
cnts s“^ MeV - ^ cm"^ compared with 0.067 cnts s“* MeV“^ cm“^ for the planar 
array. The planar array has 5 times less aperture flux background at this 
energy, which is consistent with the aperture solid angles of the two 
configurations (0.08 sr - coaxials, 0.015 sr - planars; see Table 1). Below 
100 keV, the shape of the coaxial and planar aperture flux backgrounds in 
Figures 1 and 2 are quite different from one another. The planar background 
has a maximum at 35 keV due to the peak in the incident gamma-ray spectrum at 
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this energy (Figure 3), whereas the coaxial background has a maximum at a 85 
keV due to the Ge deadlayer on top of the detectors. 

An interesting question that can be answered with the present calculation 
results Is: How small an aperture is required to reduce the aperture flux 

background to the level of the elastic neutron scattering background in new 
Instruments? The comparislon will be done at 40 keV and at an atmospheric 
depth of 3.5 g cm" 2 over Palestine, Texas. For an instrument with a 15.2 cm 
thick Nal shield, the elastic neutron scattering background Is 5.0 x 10" 2 cnts 
s"l MeV" 1 cm" 3 at 40 keV. The incident photon flux at this energy is 17 
photons cm" 2 sr" 1 s" 1 MeV" 1 or 5.2 x 10" 3 photons cm" 2 (square degree)" 1 s" 1 
MeV" 1 , and the detector efficiency is near unity. Therefore, we find that, 
for an instrument with an aperture solid angle of fi, detectors of thickness d, 
and no deadlayer on the detectors (i.e., either a planar detector or an n-type 
coaxial detector), the aperture flux background is equal to the elastic 
neutron scattering background at 40 keV if 

^ « 10 (square degrees) cm" 1 . (5) 

If the aperture size is being chosen to minimize the instrument background, 
then equation (5) gives an approximate lower limit that need be considered for 
n; reducing n much more gives diminishing returns since the elastic neutron 
scattering then dominates the background. 

For a planar detector of 1 cm thickness or a segmented coaxial detector 
(see Section 4.3) with top segment of 1 cm thickness, and assuming a square 
field of view, equation (5) recommends a field of view of size approximately 
3° x 3°. For a 5 cm, unsegmented detector, the field of view need not be much 
smaller than 7° x 7° if background reduction is the goal. 
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4.3 Detector Segmentation 

The understanding of the nature and relative intensities of the various 
background components gained with the present background study, has led to a 
new detector concept that promises to improve the sensitivity of future gamma- 
ray spectrometers. In this section, we describe the detector and calculate 
its performance as a function of various detector and instrument parameters. 

At low energies, the background in future instruments with narrow fields 
of view will be largely detector-volume dependent while the signal from 
incident low-energy gamma rays on the top surface of the detector is area 
dependent. The one background component that Is not volume dependent is the 
aperture flux, but for narrow field-of-view instruments {unlike LEGS) this 
component is small. This is Illustrated in Figure 10a which shows the 
calculated background components for an instrument with a 3° x 3° field of 
view, coaxial n-type germanium detectors of size 7 cm diameter by 7 cm length, 
a 15,2 cm thick Nal shield with threshold set at 20 keV, and 500 g per 
detector of passive A1 inside the shield. As mentioned in Section 4.1, n-type 
Ge detectors have the thin contact on the outside surface, thereby eliminating 
the thick deadlayer on the top of the detector and allowing low-energy photons 
(10 to 100 keV) to be detected with essentially no attenuation. Since the 
background below 100 keV in Figure 10a is dominated by the volume-dependent 
elastic neutron scattering component, the desire is for a thin large-area 
detector in this energy range. At high energies, however, both the background 
and the signal are volume dependent, so that the instrument sensitivity 
increases with increasing detector volume. With the LEGS instrument, these 
two opposing conditions were met by having two detector arrays - the planar 
array for low-energy observations and the coaxial array for high-energy 
observations. 
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The new idea is to have a single large coaxial detector with Its outer 
side contact segmented into several horizontal rings. Although the detector 
is maintained as a single device, the use of the signals from the segmented 
contact give It characteristics similar to those of a stack of planar 
detectors. The high-resolution signal Is still obtained from the unsegmented 
Inner contact, while the information from the segmented outer contact is 
stored as tags for each event, to be used during data analysis. For the 
present discussion of the low-energy response, the relevant segment is the top 
one. If only those events which have signal in the top segment and none in 
the lower segments are used at low energies, the background is reduced by 
approximately the ratio of the top segment volume to the total detector 
volume. At the same time, the efficiency for detecting low-energy gamma rays 
Incident through the aperture is essentially unchanged. The sensitivity in 
this energy range is therefore improved, as will be discussed more 
quantitatively below. The background reduction for a 1 cm thick top segment 
is illustrated in Figure 10b. The background per unit volume of Ge for this 
top segment mode ip, seen to be roughly the same Intensity as the volume- 
normaliy.ed all -events background in panel (a). Therefore, since the volume of 
the top segment is a factor of v 7 smaller than the volume of the whole 
detector, the background per detector is significantly reduced. 

At higher energies, the segmented detector can be used in a different 
mode to reduce the background. The idea here is that between 150 keV and 1 
MeV the background is dominated by the localized 0“ decays (Figure 10a), which 
are predominately decays to the ground state of the daughter nuclide. As 
discussed in Section 3, the signal for each event is caused by the energy loss 
of the 0“ electron, and is therefore produced in a very small region (< 1 mm) 
of the detector due to the short range of electrons at these energies. The 
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elastic neutron scattering background, which dominates below 100 keV and is 
important up to 300 keV, is also a localized interaction with ionization limited to 
an extremely small region for each event. In contrast, the incident gamma 
rays in the 150 keV to 1 MeV range interact predominantly via Compton 
scatterings, which produce ionization at more than one site in the detector. 

Hence, if only those events which have signal in more than one segment are 
accepted, the background is substantially reduced while the gamma-ray 
detection efficiency remains high. This concept is illustrated in Figure 11. 

The background for multiple-segment events is shown in Figure 10c. The 
plot is for detectors with 7 segments, each 1 cm thick, assumed to have no 
dead region between segment boundaries and to have very low segment thresholds 
(any energy deposition causes the segment to be included in the coincidence). 

With the multiple-segment requirement, the elastic neutron scattering 
component is eliminated from the background, and the aperture flux cuts off 
sharply below 100 keV since events with lower energies are almost all top- 
segment-only events. The localized 3" decays are reduced by more than an 
order of magnitude, but are not entirely eliminated because some of the 3" 
electrons cross segment boundaries. To determine the magnitude of this 
effect, a Monte Carlo program was written that propagates electrons in a 
segment. It was assumed that all electrons at a given energy have the same 
range, which was taken to be their extrapolated range (see footnote in Section 
3.3). Since the extrapolated range is at the upper end of the range 
distribution for actual electrons, this assumption results in an overestimate 
of the number of localized 3" decays that contribute to the multiple-segment 
background. The nonlocal ized 3“ decays are an important component of the 
multiple-segment background, although their Intensity is less than in the all- 
events mode due to the fact that many of the prompt gamma-rays are absorbed in 
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the same segment as the &“-decay event. The shield leakage component Is also 
less Intense and cuts off more steeply toward lower energies than In the all- 
events mode, but Is the dominant component at almost all energies. The total 
background In the multiple-segment mode Is seen to be reduced from the all- 
events mode by more than an order of magnitude between 150 and 500 keV and a 
factor of five at 1 MeV. 

The idea of segmenting the contact on Ge detectors for medical Imaging 
purposes [47,48] and for some limited background rejection applications [49- 
51] has been around for many years. Also, the concepts of using thin, large- 
area detectors to minimize background In the hard X-ray energy range and 
Compton telescopes to minimize background In the MeV energy range are not 
new. What is new In the detector concepts presented recently by Gehrels et 
al . [13], Roth et al. [52], and Varnell et al. [53], and analyzed in detail 
here, is the idea of using segmentation to reduce both the low-energy and 
medium-energy backgrounds in a single Ge detector. 

The relevant quantity for comparing the performance of a segmented 
detector with that of an unsegmented detector Is not the background level, but 
rather the instrument sensitivity for detection of lines in a source 
spectrum. The sensitivity Is the minimum flux in a spectral line that can be 
detected at a given significance level in a given amount of time, so that 
lowering the sensitivity represents an improvement. The sensitivity, S, is 
given by 


S = 2 k ( n V ae B f t) 1/2 
exp(-ux) n A e f t G 


photons cm 


-2 



( 6 ) 


where k is the significance level of the line search (# of cr), n Is the number 
of detectors of volume V {cm 3 } and area A {cm 2 }, AE is the energy Interval 
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(MeV) searched, B Is the background (cnts s"* MeV"* cm" 2 ) in that interval, f 
is the livetime fraction during the observation, t is the total observing time 
(s) assumed to be equally divided into source and background observations, u 
is the attenuation factor (cm 2 g - *) in air, x is the atmospheric depth (g cm -2 ) 
along the line of sight of the telescope pointing direction, e is the full- 
energy-peak gamma-ray detection efficiency, and G is the fraction of the line 
flux that Is contained in the Interval AE given the Instrumental resolution. 

It was assumed in the derivation of the equation that the observation is 
background limited and that the number of background counts in AE is much 
greater than 1 so that Gaussian statistics apply. In the analysis below, we 
use the following values: k=3, n a 7, ae=4x 10” 2 MeV which is close to the 
optimum energy interval for a narrow-line search [3] given the v 3 keV energy 
resolution expected for the new detectors, G=0.88 based on this choice of A! E, 
x 3 3.5 g cm” 2 , f=0.9, and t=2.88xl0 4 s (8 hours). 

An interesting aspect of both the top-segment and multiple-segment modes 
is that not all events are used when these modes are employed. A small 
additional improvement in sensitivity can therefore be obtained by also 
including in the analysis the leftover events, albeit with their higher 
background level and lower gamma-ray efficiencies. For example, in the case 
of the multiple-segment mode, single-segment events can also be analyzed and 
the two sensitivities combined in quadrature. Denoting the efficiency, 
background, and volume for the primary mode as b^, and V^, and for the 
remaining events as eg, Bg, and Vg, and using the above values for the various 
parameters in equation (6), the sensitivity becomes 

-3 e 2 e 2 

c _ 1.01 X 10 r 1 x 2 v-1/2 

5 " expT r u T 375TT + Y^> 


( 7 ) 
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For the all -events mode, there are no leftover events and the second term does 
not apply U 2 « o). 

Figure 12 shows the efficiencies for the three analysis modes obtained 
with the Monte Carlo code. Using this date, the background data in Figure 10, 
and equation (7), we have calculated the narrow-line sensitivity for the three 
segment coincidence modes assuming an instrument configuration with 7 
detectors of size 7 cm diameter and 7 cm height, each with its cathode divided 
into 7 segments, and surrounded by a 15.2 cm thick Nal shield. The results 
are shown in Figure 13. Since these sensitivity curves are based on the 
calculated continuum backgrounds not including the background lines, they do 
not apply at the energies of the strong background lines, such as 23, 67, 140, 
198, 511, 844, 1369, and 1461 keV. The solid line gives the all-events 
sensitivity which would be obtained with unsegmented detectors, and is 
approximately a factor of 4 to 6 better than existing instruments such as 
LEGS. The dashed lines show the improvement that can be obtained with a 
segmented detector using the top-segment and multiple-segment modes, and the 
upper panel gives the ratios of the all-events sensitivity to the 
sensitivities for the segmented modes. The conclusion is that the sensitivity 
can be improved by approximately a factor of 2 between 20 and 700 keV by using 
a segmented detector. The significance of a factor of 2 gain in sensitivity 
is best appreciated by considering that the number of unsegmented detectors in 
an instrument would have to be increased by a factor of 4 to achieve a similar 
gain. In the next several paragraphs, we explore the dependence of the 
sensitivity improvement on the number of segments per detector, the detector 
size, and the shield parameters. 

Figure 14 shows the sensitivity improvement factor at three different 
energies for segmented detectors (relative to unsegmented detectors) as a 
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function of the number of segments equally dividing the outer side contact of 
the detector. Panel (a) Is for 30 keV where the top-segment mode applies, and 
therefore pertains only to the thickness of the top segment. The Improvement 
at this energy Is a strong function of segment thickness and only begins 
leveling off for extremely thin top segments In the 1 to 2 mm range. The 
sensitivity gain for these thicknesses Is more than a factor of 3. Of course, 
at slightly higher energies, the incident gamma rays penetrate deeper Into the 
detector and the improvement curves level off at larger top-segment 
thicknesses. For Instance, already at 40 keV the curve levels off closer to a 
1 cm segment thickness, as was predicted by equation (5). 

Panels (b) and (c) of Figure 14 show the improvements at 0.2 and 1 MeV 
for the multiple-segment mode. In these cases, only a small gain Is achieved 
by having more than approximately 10 segments per detector. Toward smaller 
numbers of segments, the improvement at 0.2 MeV falls off steeply below 4 
segments. An Interesting special case shown by + 's In the figure Is a 2- 
segment detector with unequal segment thicknesses, 1 cm on top and 6 cm on the 
bottom. At 0.2 MeV, this detector performs approximately as well as the 7- 
segment detector, although at 1 MeV its relative performance is down. The 
reason it does so much better at 0.2 MeV than the detector with 2 equal 
segments Is that, for top-incident gamma rays at this energy, a very common 
event Is one that interacts via a Compton scattering near the top of the 
detector with the scattered photon traveling on the order of a centimeter 
before being photoabsorbed. These are Included as valid multiple-segment 
events for the 1 cm/6 cm detector but not for the equal-segment detector. At 
1 MeV, the scatterings occur typically deeper In the detector, and the 
relative performance of the 1 cm/6 cm detector is not as good. Since any 
multi segment detector with a 1 cm top segment gives the same segment 
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information as the 1 cm/6 cm detector, plus additional information, It is 
possible to devise analysis schemes for a multisegment detector that are more 
complicated than the multiple-segment mode and that take advantage of the 
specific interaction mechanisms at various energies. For Instance, at 0.2 MeV 
the multiple-segment events that Include the top segment can be analyzed 
separately from those that do not, and the sensitivities combined as in 
equation (7). The Improvement factor for a 7-segment detector Increases from 
1.9 to 2.2 using this technique at 0.2 MeV, as shown by the open circle In the 
figure. 

There are undoubtedly other segment coincidence conditions that give even 
larger sensitivity gains in specific energy ranges. The advantage of 
recording all events during the observations, regardless of segment 
coincidence, is that the optimum mode in each energy range can be chosen 
during the data analysis to give the best sensitivity. The larger the number 
of segments, the more possibilities are available, so that a trade-off exists 
between complexity in data analysis and Instrumentation versus possible 
sensitivity gains. Figure 14 indicates that diminishing returns occur for 
detectors with more than 10 segments. The figure also shows the benefits at 
low energies of having as thin a top segment as possible. 

In Figure 15, we show the effect of detector size on the sensitivity 
improvement obtainable with a segmented detector. A constant segment 
thickness of 1 cm was assumed In the calculation, so that the number of 
segments per detector increases from 5 for the 5 cm detector to 8 for the 8 cm 
detector. At 1 MeV, the sensitivity Improvement is fairly flat, Indicating on 
the one hand that segmentation is a useful technique even if one is using only 
5 cm detectors, but on the other hand that the improvement factor will not 
increase much as one obtains larger detectors. At 0.2 MeV, the difference 
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between small and larye detectors Is larger, but the general conclusions are 
the same. 

In Figure 16 the effect on the sensitivity Improvement of changing the 
shield thickness and threshold Is shown. The shield thickness Is seen to be a 
critical parameter for detector segmentation. This Is not surprising since 
the background in the multiple-segment mode Is dominated by the shield leakage 
component (see Figure 10c). For a 10 cm thick Nal shield, detector 
segmentation is probably not worthwhile for achieving sensitivity improvement 
in the multiple-segment mere. Segmentation becomes useful for shields thicker 
than ~ 12 cm, and is a powerful technique for shields in the 15 to 20 cm 
thickness range. The shield threshold level Is also an Important parameter, 
as can be seen in panels (c) and (d). The benefit of segmented detectors Is 
significantly decreased for threshold levels over v 100 keV. 

5. Conclusions 

The necessary analytical tools and data on interaction cross sections and 
atmospheric radiation fluxes are now available for performing detailed 
calculations of the background observed by gamma-ray spectrometers flown in 
the upper atmosphere. For spectrometers with Ge detectors and thick active 
shields, the components of the continuum background are: 1) elastic neutron 

scattering which is due to atmospheric neutrons that penetrate the shield and 
scatter elastically on the Ge nuclei In the detector, and is important at low 
energies ($ 100 keV); 2) aperture flux which is due to atmospheric and cosmic 
gamma rays that enter the aperture of the instrument and which can be 
important at low energies (< 100 keV) depending on the instrument field of 
view; 3) 3“ decays which are due to protons and neutrons produced in the 
atmosphere and in the shield that interact with the Ge nuclei to produce 3“- 
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unstable nuclides, and are Important between 100 keV and 1 MeV; and 4) shield 
leakage which Is due to the small fraction of atmospheric gamma rays that 
penetrate the shield without depositing more energy In It than the shield 
threshold, and Is the dominant component at high energies {> 1 MeV). 

The detailed stud|y of the background components presented In this paper 
has led to a better understanding of their nature and to several specific 
recommended techniques for reducing the background In future instruments. In 
present Instruments such as the Goddard LEGS Instrument, the aperture flux and 
shield leakage backy”ounds are Increased by factors of > 2 by passive material 
Inside the shield (such as the Ge deadlayer on the detectors, the detector 
coldflnger, the detector housings and cryostat, passive collimators in the 
aperture, and housings on the shield) and by relatively high (> 100 keV) 
thresholds in the shield. Concerning the passive material, recommendations 
include using n-type Ge detectors with the thin outer deadlayer, building 
detector housings and cryostats with the minimum possible low-Z material, 
replacing the standard Cu coldflnger with an A'l coldfinger, and using active 
collimators when possible. To obtain as low a shield threshold as possible, 
Nal is the scintillator of choice and light collection uniformity from the 
scintillator is a critical parameter. 

A new type of detector made from n-type Ge with its outer contact 
segmented into horizontal rings can be used to significantly reduce the 
background in future spectrometers. By using different segment coincidence 
modes in different energy ranges, the instrument's sensitivity to spectral 
lines can be improved at both low and medium energies. At low energies, the 
dominant background component in future narrow fiel d-of-view instruments will 
be elastic neutron scattering, which occurs uniformly throughout the 
detector. Incident low-energy gamma rays, on the other hand, interact 



predominantly tn the top region of the detector, so that by analyzing only 
top-segment events at low energies the signal to background ratio is 
Improved. For instance, for detectors of size 7 cm diameter by 7 cm height 
with a 1 cm thick top segment, inside a 15.2 cm thick Nal shield and with a 3° 
x 3° field of view, the sensitivity is improved by a factor of ~ 2 between 20 
and 100 keV. 

At medium energies (100 keV to 1 MeV), background for the above 
Instrument configuration is dominated by the 3"-decay component, which our 
calculations show is made up almost entirely of decays to the ground state of 
the daughter nuclides. Since no prompt gamma rays are emitted, and since the 
range of the 3" electron is small at these energies, ionization in the 
detector is confined to small localized regions. On the other hand, incident 
gamma rays at these energies interact predominantly via Compton scattering 
that deposit ionization in more than one region of the detector. Therefore 
the signal to background ratio is increased by accepting only multiple-segment 
events. For a 7-segment detector, this technique gives a sensitivity 
improvement of a factor of 1.5 to 2 between 150 keV and 1 MeV. The 
sensitivity improvement that can be obtained in this multiple-segment mode 
depends critically on shield thickness and threshold. For example, concerning 
shield thickness, the Improvement factor at 0.2 MeV is only 1.3 for a 12 cm 
thick Nal shield, but increases to > 2 for > 15 cm. 
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Appendix 

For each interaction contributing to the P"-decay background, three 
different nuclides are involved. The first, which will be labeled "1", is the 
Ge isotope that the incident neutron or proton interacts with. The second, 
labeled "j", is the product nuclide created in the interaction, which in turn 
beta decays to the daughter nuclide, labeled "k". The daughter nuclide may be 
produced in an excited state that, in most cases, decays promptly to the 
ground state emitting one or more gamma rays. The equation for calculating 
the count rate per unit volume of Ge for each interaction and decay is 


10- 24 pN. “ , , 

R ijk = a 1 b j d jk E k (1 “ 2 — w — - t Q °1j J dE cnts s cm 


( 8 ) 


where a-} is the abundance fraction (atoms) of the I** 1 isotope of Ge, bj is the 
P" branching fraction of the interaction product, dj k is the fraction of the 
0” decays that go to the nuclear state of interest of the daughter nuclide, e|< 
is the probability that the prompt gamma rays produced in the decay do not 
escape the detector and cause a shield veto (e k = 1 for ground-state decays), 
t is the time since the interactions started which we take to be the time 
since the instrument ascended through the Pfotzer maximum in the proton and 
neutron fluxes (t = 6 hours 1 'n this analysis), x. is the 3”-deca,y half-life of 

J 

the interaction product, 10 -24 is the number of cm 2 per barn, p is the density 
of the target material, N A is 6.02 x 10 23 atoms mole -1 , w is the atomic weight 
of the target material, Pjj is the cross section in barns for a neutron or 
proton incident on a nucleus of the i^ isotope of Ge to produce the 
unstable product, and j is the neutron or proton flux at the detector in units 
of cm -2 s -1 MeV -1 (see Figure 4). 

The rates for all possible interactions of neutrons or protons on the 
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five stable isotopes of Ge to produce (^-unstable nuclides were either 
calculated or estimated with enough accuracy to justify neglecting them. 

There are five isotopes of Ge, two types of incident particle, over 100 
possible final states that are P" unstable [35], and in many cases more than 
one decay mode, so that well over 1000 rates had to be calculated or 
estimated. In Table 4, the relevant data for approximately 700 of these rates 
are listed. The interactions have been divided into two types: the simple 

interactions which are defined as (n , y) , (n,p), (n,2p), (n,np), (n,2n), (n,a), 
(p , y) , (p,n), (p,2n), (p»np), (p,2p), and (p,a); and the complex interactions 
{sometimes called spallation interactions) which include all others. All 
simple interactions producing ^"-unstable nuclides are listed in the table, 
with the exception of the (n,2.p) and (p,y) interactions whose cross sections 
were not available. The reason we could not find cross sections for these two 
interactions is, almost certainly, that the cross sections are extremely 
small. For instance, for Mn where the cross sections are available, the value 
for (n,2p) at 15 MeV is <_ .3 millibarn [54] and for ( p , y) at 10 MeV is .1 
millibarn [55], compared with more typical values for the other interactions 
of 10 to 100 millibarn. The decays in the table with rates greater than 
roughly 1.5 x 10" 6 s"* cm" 3 are numbered, whereas those with smaller rates are 
labeled "s" for small. This cutoff is three orders of magnitude smaller than 
the highest rate of 3.9 x 10* 3 s“* cm" 3 for the ground-state decay of 75 Ge 
produced by 74 Ge(n,y) 75 Ge. The interactions labeled as X(P")Y are for the 
two-step process where a 3"-unstable nuclide is produced by the decay of 
another unstable nuclide. 

For the complex interactions in Table 4, the rates for neutrons and 
protons incident on all five Ge isotopes to yield a given product nuclide were 
summed together and listed under that nuclide. These listings are therefore 
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typically each the sum of 10 rates. In general, cross sections were not 
available for the neutron-induced complex Interactions, so we used proton 
cross sections in these cases. This should not Introduce large errors In the 
results since the neutron and proton cross sections are similar at the high 
energies required to produce the complex Interactions, 

For each decay, Table 4 indicates whether it is localized or nonlocallzed 
(see Section 3.3), and gives the 3" half-life and branching percentage, the 
decay mode and decay forblddeness, the P~ endpoint energy, the rate, and the 
cross-section reference. All of the decay data are from the Table of Isotopes 
[35]. Only those nonlocalized decays accompanied by one prompt gamma ray are 
included in the decay mode listings (except in a few cases where a second 
gamma ray has very low energy). It was assumed that, for multiple gamma-ray 
decays, the probability of having radiation escape the detector and interact 
in the shield is high. For decays with single prompt gamma rays, the factor 
e k in equation (8) was assumed to be the probability that the gamma ray is 
totally absorbed in the detector, as discussed in Section 3.3. The value of 
this probability was determined with a Monte Carlo program, and Is shown as a 
function of gamma-ray energy in Figure 17 for two different detector sizes. 

The absorption efficiencies for the LEGS coaxial detectors were used in 
calculating the rates listed in Table 4. 

Given the rates, the next step Is to determine the 3" energy spectrum for 
each decay, and then to sum the spectra into final localized and nonlocalized 
background spectra such as those shown in Figures 1, 2 and 10. The formula 
for the 3“ energy spectrum is given by Behrens and Szybisz [36] as 

N(E) a N(W) a p W (W 0 ~W) 2 F(Z,W) C(W) 


(9) 
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where E Is the kinetic energy of the electron, W = E/mc 2 +l Is the total 
electron energy In units of me 2 ,, me 2 = 0.511 MeV Is the electron rest mass, 

N(E) and N(W) are the number of decays per unit time per unit kinetic energy 
and total energy respectively, p ° (W 2 -l)^/ 2 Is the electron momentum In units 
of me, W 0 = E 0 /mc 2 + 1 Is the endpoint total energy In units of me 2 , E 0 Is the 
endpoint kinetic energy (Table 4), F(Z,W) Is the Fermi function, Z Is the 
atomic number of the daughter nucleus, and C(W) Is a correction factor whose 
form depends on the forblddeness of the decay (Table 4) and whose numerical 
parameters are given by Behrens and Szyblsz. We obtained values for the Fermi 
function and other functions needed for calculating C(W) from the tandolt- 
Bffrnsteln data tables [56]. The spectra were normalized to have Integrals 
equal to the rates In Table 4. For those decays with a prompt gamma ray, the 
3“ spectrum was offset by the gamma-ray energy since both signals are 
collected in the detector. 

The 3" decay spectra calculated using equation (9) and the data In Table 
4 are shown in Figure 18 for all the decays that are numbered in the table. 

The figure is divided into three panels In order to minimize confusion among 
the spectra. The localized decays have short-dashed curves, and the 
nonlocalized decays have long-dashed curves. The solid line is the sum of the 

localized decays. For the localized decays, the dominant interaction Is 

7 ^Ge(n,Y) 75 Ge followed by the complex interactions producing 70 Ga and 6q Zn. 
Although all other decays fall considerably below these three, there are a 
large number of them, and their contribution makes up - 20% of the total below 

60 keV and over half of the total above 1 MeV. For the nonlocalized decays, 

74 Ge(n,Y) 75 Ge is by far the dominant interaction. 



48 


References 


[1] W, Paciesas, R. Baker, D. Boclet, S. Brown, I. Cline, H. Costlow, Ph. 
Durouchoux, C. Ehrmann, N, Gehrels, J. M. Hameury, R. Haymes, 8. 
Teegarden, and J. Tueller, Nucl. Instr. and Meth. 216 (1983) 261. 

[23 E, A. Womack and J. W, Overbeck, J. Geophys. Res. (1970) 1811. 

[3] A. S. Jacobson, R. J. Bishop, G. W. Culp, L. Jung, W. A. Mahoney, and J. 
B. Willett, Nucl. Instr. and Meth. 127 (1975) 115. 

[4] W. L. Imhof, G. H. Nakano, and J. B. Reagan, J. Geophys. Res. jU (1976) 
2835. 

[5] M. Leventhal, C. MacCallum, and A. Watts, Astrophys. J. 216 (1977) 491. 

[6] W. A. Mahoney, J. C. Ling, J. B. Willett, and A. S. Jacobson, In Gamma 
Ray Spectroscopy In Astrophysics , eds. T. L. Cline and R. Ramaty, NASA 
TM-79619 (1978) 462. 

[7] F. Albernhe, D. Boclet, J. P. Chabaud, J. Clalsse, Ph. Durouchoux, P. 
Frabel, E, Olivier, J. M. da Costa, P. Pagnler, R. Rocchia, and G. 
Vedrenne, Nucl. Instr. and Meth. 155 (1978) 171. 

[8] L. E. Peterson, R. M. Pelllng, and J. L. Matteson, Space Sci. Rev. 13_ 
(1972) 320. 

[9] J. C. Ling, Thesis, Univ. of California, San Diego (1974). 

[10] J. L. Matteson, P. L. Nolan, W. S. Paciesas, and R. M. Pelling, Space 
Sci. Instr. 3 (1977) 491. 


49 


[11] G. R. Ricker, J. V, Vallerga, and D. R. Wood, Nucl. Instr. and Meth. 213 
(1983) 133. 

[12] J. C. Ung, W, A. Mahoney, J, B. Willett, and A. S. Jacobson, J. 

Geophys. Res. 82_ (1977) 1463. 

[13] N. Gehrels, T. L. Cline, B. J. Teegarden, J. Tueller, M. Leventhal, C. 

J. MacCallum, P. V. Hewka, and P. Ryge, IEEE Trans. Nucl. Scl. NS-31 , 

No. 1 (1984) 307. 

[14] R. L. Klnzer, private communication (1983). 

[15] R. L. Klnzer, W. N. Johnson, and J. D. Kurfess, Astrophys. J. 222 (1978) 

370, 

[16] V, Schonfelder, F. Graml, and F.-P. Pennlngsfeld, Astrophys. J. 240 
(1980) 350. 

[17] J. A. Lockwood, L. Hsleh, L. Frlllng, C. Chen, and D. Swartz, J. 

Geophys. Res. 8£ (1979) 1402. 

[18] J. C. Ling, J. Geophys. Res. 80 (1975) 3241. 

[19] J. M. Ryan, B. Dayton, S. H. Moon, R. B. Wilson, A. D, Zych, and R. S. 

White, J. Geophys. Res. j82 (1977) 3593. 

[20] R. S: White, B. Dayton, S. H. Moon, J. M. Ryan, R. B. Wilson, and A. D. 
Zych, Astrophys. J. 218 (1977) 920. 

[21] V. Schffnfelder, U. Graser, and J. Daugherty, Astrophys. J. 217 (1977) 
306. 

[22] R. D. Evans, The Atomic Nucleus, McGraw-Hill, New York (1955) 625. 


50 


[23] K. F. Koral and A. J. Cohen, NASA TN D-2909 (1965). 

[24] R. L. Kelly et a!.. Rev. Mod. Phys. .52, No. 2, Part II (1980) SI. 

[25] T. W. Armstrong, K. C. Chandler, and J. Barish, J. Geophys. Res. 78_ 
(1973) 2715. 

[26] B. Rossi, High-Energy Particles , Prentice-Hall, Englewood Cliffs, NO 
(1952) 486. 

[27] F. B. McDonald and W. R. Webber, Phys. Rev. 115 (1959) 194. 

[28] W. R. Webber and J. A. Lezniak, Astrophys. and Space Sci. 3CM1974) 361. 

[29] W. E. Alley and R. M. Lessler, Lawrence Livermore Lab. Report No. UCRL- 
5048s, rev. 1 (1972). 

[30] Handbook on Nuclear Activation Cross-Sections , Inti. Atomic Energy 
Agency Tech. Report Series No. 156, Vienna (1974). 

[31] D. I. Garber and R. R. Kinsey, Brookhaven National Lab. Report No. 325, 
3rd ed., vol. II (1976). 

[32] R. J. Howerton, D. E. Dye, and S. T. Perkins, Lawrence Livermore Lab. 
Report No. UCRL-50400, vol. 18, rev. 1 (1982). 

[33] R. Sllberberg and C. H. Tsao, Astrophys. 0. Suppl.^25 (1973) 315. 

[34] , Astrophys. J. Suppl. £5_ (1973) 335. 

[35] Table of Isotopes , eds. C. M. Lederer and V. S. Shirley, 7th ed. , Wiley, 
New York (1978). 



[36] H. Behrens and L. Szybisz, Shapes of Beta Spectra, Phys. Data 6-1 
(1976). 

[37] G. F. Knoll, Radiation Detection and Measurement , Wiley, New York (1979) 
571. 

[38] D. I. Garber, L. G. Strffmberg, M. D, Goldberg, D. E, Cullen, and V. M. 
May, Brookhaven National Lab. Report No. 400, 3rd. ed., vol. I and II 
(1970). 

[39] 0. W. Corbett, Electron Radia tion Damage In Semiconductors and Metals , 
Academic Press, New York (1966) 135. 

[40] J. Llndhard, V. Nielsen, M. Scharff, and P. V. Thomsen, Danske Vid. 
Selskab, Mat.-fys. Medd. 33, No. 10 (1963). 

[41] M. T. Robinson, In Nuclear Fusion Reactors , eds. J. L. Hall and J. H. C. 
Maple, British Nucl . Energy Soc., London (1970) 364. 

[42] R. R. Daniel and S. A. Stephens, Rev. Geophys. and Space Phys. 12 ^ (1974) 
233. 

[43] J. M. Lavlgne, M. Niel, G. Vedrenne, B. Arglnler, E. Bonfand, B, 

Parlier, and K. R. Rao, Astrophys. J. 261 (1982) 720. 

[44] G. Fulks and P. Meyer, J. Geophys. _40 (1974) 751. 

[45] J. F. danni, Air Force Weapons Lab. Tech. Report No. AFWL-TR-65-150, 
Kirtland Air Force Base, New Mexico (1966). 

[46] T. W. Raudorf, R. C. Trammell, and L. S. Darken, Jr., IEEE Trans. Nucl. 
Sc 1 . NS-26, No. 1 (1979) 297. 


[47] R. P. Parker, E. M. Gunnersen, J. L. Wankling, and R. Ellis, Medical 
Radioisotope Scintigraphy _1_, IAEA (1969) 71. 

[48] M. G. Strauss and I. S. Sherman, IEEE Trans. Nucl. Sci. NS-19, No. 3 
(1972) 219. 

[49] C. R. Gruhn, J. V. Kane, W. H. Kelley, T» Kuo, and G. Berzins, Nucl. 
Instr. and Meth. 54_ (1967) 268. 

[50] J. M. Palms, R. E. Wood, and 0. H. Puckett, IEEE Trans. Nucl. Sci. NS^ 
15, No. 3 (1968) 397. 

[51] A. R. Sayres and J. A. Baicker, IEEE Trans. Nucl. Sci. NS-15 , No. 3 
(1968) 393. 

[52] J. Roth, J. H. Primbsch, and R, P. Lin, IEEE Trans. Nucl. Sci. NS-31 , 
No. 1. (1984) 367. 

[53] L. S. Varnell, 0. C. Ling, W. A. Mahoney, A. S. Jacobson, R. H. Pehl, F 

S. Gouldlng, D. A. Landis, P. N. Luke, and N. W. Madden, IEEE Trans. 

Nucl. Sci. NS-31 , No. 1 (1984) 300. 

[54] E. T. Bramlitt and R. W. Fink, Phys. Rev. 131 (1963) 2649. 

[55] S. E. Woosley, W. A. Fowler, J. A. Holmes, and B. A. Zimmerman, 

unpublished preprint OAP-422, Calif. Inst. Tech. (1975). 

[56] H. Behrens and J. Jffnecke, in Landolt-Bffrnstein Numerical Data and 
Functional Relationships in Science and Technology , New Series, Group I 


Vol . 4, Springer, New York (1969). 


LEGS Observations Used for Background Calculations 




f’S- 

C OJ 
QMU 01 
12 M U 

•r* C7l 

P E OJ 
C p TO 
*r- U ' — ' 

O 4 - 

a. 


u. o 

P ii* 

is u — 
at at jtcvj 
E£*J! 
5 aap 
fc to at o 
■MOO 


3 C 
i — O <u 
(OLE 
> o) a* *- » 

i- M t. H- 
at o 3 sa 
P ra to — « 
c ca its 

t-t QJ 

t- e: 

QJ O 

6 U. 


X5 "O 

i— TO r— 
Hr O** 
e dj jc > 

•r* 'p- tfl 01 

E£ 0)^ 
o in u — 

= jE 


CD 

ft* r— Vi 

qj a> c 
i- e re 

3 ^ *f- 
P *3 
i. TO <0 
QJ >r- U 
ELr- QJ 

< O P 
</> to 


>-> 5 - 



4-> 


o 





3 


p 





O 


o 





J= 


no 





O) 


4- 





3 






o 

o 


re 




O 

t- 


10 




H 

J= 






+-> 


re 





cn 


XT 





e 


u 





•r- 


3 





4-> 


E: 





O 






CM 

03 


tfl 




1 

i~ 


re 




o 

0) 






vH 

4-J 


c 


• 


X 

C 


at 


p 


LO 

♦r* 


QJ 


c 


* 



jo 


QJ 


*H 

(0 




E 



C 


QJ 


QJ 



o 


> 


t- 



p 


re 


3 



o 


JO 


to 

■ 


x: 



• 

re 

•to 


CL 


>* 

P 

aj 

r~ 




S 

c 

E 

at 

J- 

o$ 


& 

i 

TO 

*r* 

JO 

ft* 

s- 


P 

aj 

C 

to 

c 

QJ 


sz 

s_ 

3 


ft* 

c 


a i 

3 

O 

p 

pM 

QJ 


«p» 

to 

i_ 

o 

Gu 

1 


r— 

re 

D1 



J= 


4- 

at 

J > t 

c 


CD 



E 

U 

o 


•r- 


C 


re 

*r— 


J= 

* 

•r- 

TO 

sz 

P 



s. 


c 


u 


t- 

tfl 

to 

3 

cn 

QJ 


O 


T3 

o 

c 

(0 


Ll. 

<M 

f— 

s_ 

•r— 




| 

o 

D) 

J. 

" 

VI 


O 

jC 


3 

cn 

ftJ 

• 

*H 

to 

U 

•a 

,3 

X 

to 

X 

QJ 

re 


I— 

o 

St- 

CO 

U 

x> 

aj 

Q. 

CO h- 

o 

* 

J=S 


cn 

& 

p 

l»* 

P 

o> 

c 


• n 

a 



c 

re 

E 

P QJ 

a> 

to 

TO 

•f— 

t. 

o 

a. c 

p 

•r* 

i — ■ 

j. 


P 

qj .r- 

QJ 


aj 

3 

QJ 

P 

00 P 

•a 

at 

*r* 

TO 

i — 

o 

t0 


r-* 

JO 


CD 

J3 

O QJ 

4- 

O) 

V) 

aj 

C 


00 r- 

o 

C 


D) 

re 

i- 

CTi 3 


3 

i — 

c 


O 

i— t Q_ 

Q. 


rd 

re 

cn 

4- 

O 

■a 

3 

u 

e 



P 

•r* 

P 


•r* 

?- 



r— 

O 

sz 

P 

QJ 


E 

O 

re 

p 

C 



o 

(0 


£X 

•r* 



s_ 


aj 

QJ 

O 

o 


u_ 


sz 

a 

a. 

un 

5- 



h- 





p ■* 


ft « 

(> ■ * 

p**. 



rd 


JO 

u 

TO 

QJ 


54 


Table 2 

Passive Material Inside LEGS Shield 


Item 

Material 

Mass 

(g) 

Detector Configuration 

Cryostat 

A1 

695 

Coaxial 


Stainless 

930 



A1 

700 

Planar 

Detector Housings 

A1 

380 

Coaxial 

A1 

150 

Planar 

Ge Deadlayer 

Ge 

370 

Coaxial 


Ge 

125 

Planar 

Col dfinger 

CU 

55 

Both 

Shield Housing 

A1 

2000-34009) 

Both 

Passive Collimator 

Fe 

4800 

Planar 


a ^Range of values represents uncertainty in thickness of housing walls. 
Minimum thickness is estimated at 0.16 cm and maximum thickness at 0,32 cm. 



Sources of Increase In Aperture Flux Background 
for Coaxial Detectors 


Source 

Mass 

(g) 

Background 3 ) 
Increase 
(%) 

Single Scattering In Passive Materials 



Ge Deadlayer 

370 

14.9 

A1 Detector Housing 

380 

5.8 

A1 Cryostat 

695 

6.1 

Steel Baseplate . 
Al Shield Housing 5 ' 

930 

3.8 

3400 

8.0 

Multiple Scattering In Passive Materials 


13.2 

Scattering In Shield, < 200 keV 


18.8 

UE < 100 keV) 


(6.8) 

Total, 200 keV Threshold 


70.6 


a ) In 0.1 to 1.0 MeV energy range. 

b) Thick (0.32 cm) shield housing case. 


Interactions Contributing to 8“ Decay Background 
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Interaction Event 3" Branching Decay Decay Energy Rate 

Number or Product Type Half-Life Percentage Mode Forbiddeness (MeV) (Decays s'^cnfH References 
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Interaction Event B" Branching Decay Decay Energy Rate 
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Figure Captions 


Figure 1 - Observed and calculated background In the LEGS coaxial array. The 
count rate spectrum was divided by the volume of active Ge in the detectors 
(230 cm 3 ) to give units of cnts s" 1 MeV" 1 cm" 3 . The data are from flight III 
(Palestine, Texas; 5 g cm" 2 ), and the calculations are for the continuum 
background underlying the background lines. The elastic neutron scattering, 
aperture flux, 3“ decays localized and nonlocal Ized, and shield leakage 
components of the calculation are shown separately. 


Figure 2 - Same as Figure 1, but for the LEGS planar array flight V (Palestine, 
Texas; 3,5 g cm" 2 ). The Ge active volume for the planar array Is 57 cur*. 


Figure 3 - Measurements of the total downward gamma-ray flux at 5 g cm" 2 over 
Palestine, Texas. The solid line fits the data of Kinzer [14], Ki nzer et al. 
[15], and Schffnfelder et al. [16], and is the spectrum used ir/ this paper. 

The three points from Lockwood et al. [17] have been multiplied by a factor of 
1.15 [18] to correct for a depth of 3.5 g cm" 2 . The points from Ryan et al . 
[19] and White et al, [20] were measured at a zenith angle of 15°, but no 
correction was required since the zenith angle flux distribution Is 
approximately flat between 0° and 70° [18,21]. 


Figure 4 - Neutron and proton spectra inside the LEGS shield during flights 
III and V. The residual atmosphere above the Instrument plus the shield are 
equivalent to ^ 35 g cm" 2 of air. The ordinate Is the omnidirectional flux, 
which is the total flux Incident on a unit sphere. The solid lines show the 
spectra used In this paper. 
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Figure 5 - The elastic neutron scattering cross section as a function of 
Incident neutron energy for scatterings which give a Ge atom recoil energy of 
.04 MeV and .17 MeV. Also shown for aid In defining the shape of the curve at 
E n > 10 MeV, are cross sections for Cu. The solid lines are the power-law 
approximations of the data used In this paper. The data are from Garber et 
aK [38]. 

Figure 6 - The fraction of the recoil energy lost to signal-producing Ionizing 
collisions for elastic neutron scattering In Ge. The curve was calculated 
with equation (4). 

Figure 7 - The ratio of total aperture flux background to that produced by 
unscattered photons for the LEGS coaxial and planar detector arrays. 
"Unscattered photons" means those that interact only in the detector, and not 
In the passive material or in the shield. The results were calculated using 
the Monte Carlo code to simulate the Interaction of Incident photons with the 
Instrument. The error bars represent statistical uncertainties In the Monte 
Carlo results. For the coaxial array, the range in possible values due to 
uncertainties in the shield housing thickness and shield threshold are shown 
by the two curves. For the planar array, this range is smaller than the 
statistical uncertainties in the simulation and Is therefore not shown. 

Figure 8 - Total shield leakage backgrounds, including scatterings in passive 
material near the detectors and In the shield (AE < shield threshold), 
compared with the unscattered component for the LEGS coaxial array. 



Figure 9 - The effect on the shield leakage background of raising the 
threshold level in the shield. The ratio of the shield leakage background for 
a given shield threshold to the background for a 0 keV threshold is shown as a 
function of the threshold for three different energy ranges. The error bars 
represent typical absolute uncertainties in the points; the relative 
uncertainties between points on a given curve are smaller. 

Figure 10 - The calculated continuum background and its components in three 
different modes for the future Instrument configuration described in the 
text. The background is given per unit volume of Ge, with relevant detector 
volumes for each mode shown. The plotted backgrounds are for events that have 
no shield coincidences and that satisfy the segment conditions for each 
mode. The segment conditions are a) all events accepted - no conditions on 
segments, b) only events with energy deposition in top segment and none in 
bottom six segments, and c) only events with energy deposition in more than 
one segment. 

Figure 11 - Typical gamma-ray and background events in a multisegment 
detector. By requiring signal in more than one segment, the background from 
localized B“ decays is eliminated while the signal from incident gamma rays is 
kept. 

Figure 12 - Ful 1-energy-peak efficiencies for a 7 cm detector with 7 1-cm 
segments for the three segment coincidence modes. 

Figure 13 - The narrow-line ( ae= 4 keV) sensitivity for three segment modes for 
the instrument described in the text, flown at 3.5 g cm" 2 over Palestine, 


AiW * V.*.. % 
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Texas. The upper panel gives the sensitivity Improvement for the top-segment 
and multiple-segment modes compared with the all-events (unsegmented detector) 
mode. The curves do not apply at energies of strong background lines. 

Figure 14 - The sensitivity Improvement factor obtained with segmented 
detectors relative to unsegmented detectors as a function of the number of 
segments (or segment thickness) at three energies. Panel (a) is for the top- 
segment mode, and panels (b) and (c) for the multiple-segment mode. The 
instrument configuration Is the same as that assumed for Figures 10 and 13, 
with the filled circle on each curve corresponding Identically to the ratio at 
that energy in Figure 13. In all cases, except the +'s In panels (b) and (c) 
equal-size segments were assumed. The + 's correspond to a 2-segment detector 
with 1 cm top-segment thickness and 6 cm bottom-segment thickness. The open 
circles In panels (b) and (c) are the improvement factors If the multiple- 
segment events In a 7-segment detector are divided into those that include the 
top-segment and those that do not. Similar points above the curve occur for 
other multisegment detectors, but are not plotted. 

Figure 15 - The sensitivity Improvement factor obtained with segmented 
detectors relative to unsegmented detectors as a function of detector diameter 
at two energies. The detectors were assumed to be cylinders with diameter 
equal to height, and divided into 1 cm segments. The shield configuration is 
the same as that assumed for Figures 10 and 13. The filled circles on the 
curves correspond identically to the ratios at those energies In Figure 13. 

The absolute sensitivities (ph cm -2 s** 1 ) for an unsegmented detector ranges 
from (a) 1.1 x 10 -4 for 5 cm diameter to 7.8 x 10“ 5 for 8 cm at 0.2 MeV, and 
(b) 2.1 x 10“ 4 for 5 cm to 9.4 x 10~ 5 for 8 cm at 1 MeV. 
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Figure 16 - The sensitivity Improvement factor obtained with segmented 
detectors relative to unsegmented detectors as a function of shield thickness 
and threshold fir two energies. Apart from the shield, the Instrument 
configuration Is the same as that assumed In Figures 10 and 13, with the 
filled circles corresponding Identically to the ratios at those energies In 
Figure 13. The absolute sensitivities (ph cnr 2 s" 1 ) for the unsegmented 
detectors range from (a) 9.9 x 10“ 5 for 10 cm thickness to 8.7 x 10”® for 20 
cm at 0.2 MeV, (b) 1.7 x 10"4 for 10 cm to 1.1 x 10-4 f 0r 20 cm at 1 MeV, (c) 
8.6 x 10"® for 0 keV threshold to 9.1 x 10”® for 180 keV at 0.2 MeV, and (d) 
1.2 x 10-4 f 0r 0 keV to 1.3 x 10-4 f or iso keV at 1 MeV. 

Figure 17 - The probability that gamma rays produced Internal to a detector 
are fully absorbed In the detector, for two different detector sizes. The 
LEGS coaxial detector Is approximately 4.6 cm diameter by 4.6 cm length, and 
the 7 cm coaxial detector Is 7 cm diameter by 7 cm length. 

Figure 18 - Components of the P"-decay background. Identification and 
relevant data for each decay spectrum are given In Table 4. The short-dashed 
curves are localized decays and the long-dashed curves are nonlocalized. The 
sum of the localized decays is shown by the solid curve In panel (a). 
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